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AERODYNAMIC PARAMETERS OF THE NAVION AIRPLANE 
EXTRACTED FROM FLIGHT DATA 


By William T. Suit 
Langley Research Center 

SUMMARY 

An iterative method, which is characterized as a maximum-likelihood minimum- 
variance technique, was used to extract the aerodynamic parameters of a Navion airplane 
from flight data. The purposes were to compare the results with parameters obtained 
from wind-tunnel tests and with results obtained by analog matching of the same data, and 
to develop techniques for application of the parameter -extraction program. 

Results from the study showed that the parameter- extraction program can produce 
aerodynamic parameters which will permit close estimation of the aircraft time histories 
used in the extraction process. The program determined an estimate of the standard 
deviations of the states and parameters. These estimates were used to indicate how well 
the calculated states fit the flight data and the confidence in the values of the estimated 
parameters. The study also showed that the values of the parameters were affected by 
the data and mathematical model used during the extraction process. Because of the lack 
of confidence in the parameters extracted by using some of the sets of data, several 
parameters were estimated by other methods. By using a combination of methods, a set 
of parameters which gave a fit to the data was obtained. 

The extracted parameters agreed reasonably well with the values obtained by analog 
matching the same data, with the exception of the change in normal-force coefficient with 
angle of attack (CzqJ. The agreement with wind-tunnel parameters was not as good for 
the variations of pitching-moment coefficient with angle of attack (c mo! ), side-force coef- 
ficient with sideslip angle rolling- moment coefficient with sideslip angle and 

yawing-moment coefficient with sideslip angle (Cn^- However, of the parameters deter- 
mined by the program, only one had a standard deviation greater than 15 percent of the 
value of the parameter and the parameters determined gave a reasonable fit to the flight 
data. 


INTRODUCTION 

Mathematical analyses of flight dynamics and handling qualities of an aircraft are 
required for determining the suitability of the aircraft for its mission. In order to make 
such analyses, it is necessary to have available the aerodynamic parameters of the air- 
craft. There are several methods of obtaining the parameters. These methods include 



those presented in various books, wind-tunnel tests, and extraction of derivatives from 
flight-test data. Of these methods, derivatives from flight tests should be the most accu- 
rate since such results are obtained with the actual aircraft in its proper environment. 
There is, therefore, a continuing interest in developing and evaluating improved methods 
of extracting derivatives from flight data. 

In a recent study (ref. 1), a comparison was made between various analytical meth- 
ods, wind-tunnel measurements with a full-scale airplane, and results obtained from 
flight-test data for a Navion airplane. In that study, an analog-matching technique (ref. 2) 
was used in extracting parameters from the flight data. Some rather large differences 
were found between the various methods, hi particular, some large differences were 
obtained between the wind-tunnel and the flight-test results. Analog matching requires 
a highly experienced operator to match flight data properly. It appeared desirable to use 
an alternate method of extracting the derivatives from the flight data. The method used 
in this study is a mathematical formulation of the logic required to select derivatives to 
best match a set of flight data. The method selected is an iterative procedure which 
selects the aerodynamic parameters to maximize a conditional maximum likelihood func- 
tion and is equivalent to determining the set of aerodynamic parameters which will maxi- 
mize the probability that the calculated state of an airplane will match the measured state 
for the same control inputs (ref. 3). The maximization process used minimizes the mea- 
surement error covariance matrix. The resulting parameter adjustment equations are of 
the same form as those obtained by use of a modified Newton- Raphson or weighted least- 
squares technique (ref. 4). The main difference is that with the maximum likelihood for- 
mulation, the weights are updated at each iteration. The program will speed up derivative 
determination, give a fit to the flight data based on mathematically minimizing a cost cri- 
terion, and determine a matrix which indicates the variances of dependency between the 
estimated derivatives. 

The primary purpose of the present paper is to use the flight data employed in the 
analysis reported in reference 1 and extract the aerodynamic parameters for comparison 
with the results presented in reference 1. A second purpose of this paper is to indicate 
the procedure used in applying the parameter estimation program to the data herein. A 
third purpose is to relate the experience gained from this investigation and to point out the 
advantages of the program used. A fourth purpose is to indicate the confidence in the 
parameters obtained. 


SYMBOLS 

Values are given in both SI and U.S. Customary Units. The measurements and cal- 
culations were made in U.S. Customary Units. The aerodynamic parameters are refer- 
enced to a system of body axes with the origin at the aircraft center of gravity, and with 
body axes orientation as shown in figure 1. 
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acceleration, m/sec 2 (ft/sec2) 


wing span, m (ft) 

wing mean geometric chord, m (ft) 

force, N (lb) 

acceleration due to gravity, m/sec2 (ft/sec2) 
moment of inertia, kg-m 2 (slug -ft 2 ) 
index 

tail incidence angle, radians or degrees 
weighting factor 
likelihood function 

distance from aircraft center of gravity to center of pressure of horizontal 
tail, m (ft) 

moment, N-m (ft -lb) 

mass, kg (slugs) 

number of data points 

change in parameter from iteration to iteration 
rate of roll, radians/sec 
rate of pitch, radians/sec 
dynamic pressure, i-pV 2 , N/m 2 (lb/ft 2 ) 
estimate of error covariance matrix 



rate of yaw, radians/sec 


wing area, m 2 (ft 2 ) 
Thrust 


-, nondimensional 


(Dynamic pressure) (Wing area) ’ 
velocity along X body axis, m/sec (ft/sec) 

aircraft total velocity, m/sec (ft/sec) 

velocity along Y body axis, m/sec (ft/sec) 

velocity along Z body axis, m/sec (ft/sec) 

angle of attack, radians 

sideslip angle, radians 

flight-path angle, radians 

control deflection, radians or degrees 

pitch angle, radians 

air density, kg/m 2 (slugs/ft 2 ) 

roll angle, radians 

rolling-moment coefficient, M x ^qSb 
pitching- moment coefficient, My/qSc 
yawing-moment coefficient, M z yqSb 
axial-force coefficient, F x /qS 
side-force coefficient, Fy/qS 
normal-force coefficient, F z /qS 
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Subscripts: 

a aileron 

b body 


c 


computed 



e 


elevator 


f 

m 

o 

r 

t 

X 

Y 

Z 


flap 

measured 

indicates coefficient at trim conditions 
rudder 

indicates state at trim conditions 
X-axis 


Y-axis 

Z-axis 


Superscript: 

T transpose 


A dot over a symbol signifies a derivative with respect to time. 

DESCRIPTION OF AIRPLANE AND FLIGHT TESTS 


Instrumentation 

The flight data used for extracting derivatives was obtained from flight tests of the 
Princeton University variable -stability Navion airplane, N91566. The physical character- 
istics of the Navion are presented in figure 2 and table I. The data recorded for this study 
included: 



Accuracy 

Response frequency 


Hz 

cps 

Normal acceleration 

±0.01 g 

2 

2 

Roll rate 

±0.044 radians/sec 



Pitch rate 

±0.024 radians/sec 



Yaw rate ! 

±0.010 radians/sec 



Angle of attack 

±0.8° 

4 

4 

Altitude 

±30.48 m (±100 ft) 

2 

2 

Indicated airspeed 

1.03 m/sec (±2.3 mph) 

2 

2 

Control surface position 

±1 percent 

2 

2 
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The data were sequenced by a commutator at a rate of 10 points per second and 
telemetered to a ground station as a frequency-modulated signal where it was recorded 
on magnetic tape. The angle of attack was corrected for upwash effects but required no 
correction for angular rates. The accelerations and angular rates needed no correction. 

Flight Data 

The flight -test data used in the present study were obtained from tests made by pef - 
sonnel of the Princeton University Aeronautical Laboratory. The data were recorded on 
magnetic tape and were processed at the Langley Research Center for parameter extrac- 
tion. The calibrations used during processing were furnished by Princeton University. 
The data processing included digitizing the data, converting the recorded signal to engi- 
neering units, and interpolating to give all the data at the same times on the tape. The 
interpolation was required because the commutated data gave each state at a different 
time. Flight -test conditions are listed in table II. The applied control disturbances 
included an elevator doublet, an aileron doublet, and a rudder pulse. 

PARAMETER-ESTIMATION PROCEDURE 


The parameter -estimation procedure used in this study is an iterative procedure 
which maximizes the conditional likelihood function L(aerodynamic parameters, weights, 
initial conditions): 


L = 


(2,) 1 /2| R |l/2 


exp 


1 

2 


N 


Y ( Xim 

1=1 


-X ic ) R-^X; 


im 



where R is the estimate of the error covariance matrix and X is the vector describing 
the state of the aircraft. Maximizing the likelihood function minimizes the difference 
between the measured and calculated aircraft motions. 

The weighting matrix R”1 can be the complete error covariance matrix, the diag- 
onal terms of the error covariance matrix, or a diagonal matrix with fixed weights on the 
diagonal, at the discretion of the investigator. If the diagonal form of the weighting matrix 
is used, the weights represent the estimated lower bound of the noise on the measured 
states. The use of the likelihood function in parameter identification is discussed in ref- 
erence 3. Maximizing the likelihood function results in a parameter updated equation 
which is given by 

AP = (M T R-%) _1 M T X m 
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where M is the matrix of sensitivities of the calculated states with respect to the 
unknown parameters (ref. 5). The matrix will be called in this paper the estimated 
parameter covariance matrix. The updated equation is determined by forming a set of 
differential equations with the changes in the unknown parameters as the variables. This 
set of simultaneous equations is then solved by least squares to give the updated equation. 
(See ref. 5.) 

The steps in the iteration procedure are outlined in figure 3. The procedure is to 
write a set of equations of motion for the aircraft under consideration. These equations 
will include a number of aerodynamic parameters. The parameters must be initially 
estimated so that the motions of the aircraft can be calculated. These calculated motions 
are then compared with the motions of the actual aircraft for identical control inputs. 

The parameter estimation program uses the differences between the measured and cal- 
culated data to calculate the updated values of the parameters. The parameters are then 
corrected and new aircraft motions are calculated. The process is repeated with updated 
parameters until the difference between calculated and measured motions are within some 
acceptable range. The complete details are given in reference 3. 

During this investigation the mean-squared error between the measured and calcu- 
lated states was displayed at each iteration. When the mean-squared error became con- 
stant for several iterations, the problem was terminated. A printout of an estimate of the 
variance of the states, the changes in unknown parameters at each iteration, the estimated 
lower bound on the standard deviations of the unknown parameters, and the determinant of 
the R matrix were obtained. Data from the printout were examined to determine the 
fit to the flight data and the confidence in the extracted parameters. If all these criteria 
indicated a fit in the order of the instrument uncertainty and if the unknown parameters 
had changes less than 5 percent from iteration to iteration, the parameters obtained were 
considered to be as good as could be determined. 

The program is set up so that the fit to the flight data can be monitored on a cathode 
ray tube (CRT) as the parameters are updated. The mean-squared error was displayed 
on a digital voltmeter. The program can be stopped at any iteration and the states used 
in the likelihood function changed or unknown parameters added to or taken from the 
mathematical model. Also, any of the parameters in the mathematical model can be con- 
sidered to be known and fixed at specified values. These characteristics of the program 
allow the operator to have a close interaction with the program and give the operator a 
very flexible tool to use. Because the operator can see any large effect of changes in the 
parameters, there is a considerable saving in time over letting the computer run unob- 
served for a fixed number of iterations and then examining a printout of the results. Also, 
the turnaround between runs is seconds rather than hours if batch processing had been 
used. The operator’s console and cathode ray tube (CRT) are shown in figure 4. 
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PRELIMINARY CONSIDERATIONS 


Several decisions must be made prior to application of the parameter -extraction 
program. These decisions are: 

(1) What parameters should be extracted? 

(2) How good an initial estimate is required to start the procedure ? 

(3) What amount of data and what sections of the flight records should be used? 

In this study, all the parameters used in the equations of motion along with the ini- 
tial conditions on the states were extracted initially. If it was found that some of the 
parameters could not be well determined, then these were fixed at an assumed value, 
dropped from the mathematical model, or the mathematical model was changed to deter- 
mine the parameter better. The experience of the investigator in this study indicates 
that initial estimates of the aircraft parameters obtained from reference 6 gave a fit good 
enough to allow the problem to start. 

An important aspect of selecting data for analysis is that of using a section of the 
flight record which contains motion of reasonable amplitude (well above the measuring 
instrument noise level). It also is advisable to include several seconds of flight record 
at trim conditions. The flight data used in the present analysis were recorded on mag- 
netic tapes. The data for the analyses were taken from the tapes at the rate of 20 data 
points per second for each state used. 

EQUATIONS OF MOTION 

The equations of motion used herein are written relative to the body axes and are 
essentially the same as the equations given in reference 7. The assumptions used when 
writing the equation are: 

(1) The moment of inertia I^z = 0. (See ref. 8.) 

(2) The quantities q, /3, r, and p are initially zero. 

(3) The equations are basically uncoupled so that a and q variations are negli- 
gible during lateral maneuvers. 

(4) The quantities 9 and <p were initially assumed to be zero. 

(5) When 9 and cp were used in the equations, they were calculated by using 
9 = q cos (p - r sin (p and (p = P + (q sin <p + r cos $)tan 9. 

(6) The velocities u, v, and w were calculated from the a, p, and V flight 
data by u = V cos a cos p, v = V sin / 3 , and w = V sin a cos p. 
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(7) In the v equation, u and w were assumed to be constant and fixed values 
were put in. 

By using these assumptions, the equations of motion used were for the longitudinal 
equations: 


ii = -qw - g sin 6 + | p ~ ^jc x>0 + C Xff (a - a t j] 
w = g cos 0 + qu + ip ^[c z ,o + Cz a { a - «t) + c Zg e (® e “ 6 e t )J 

|c m ,o+ S(“ “ a t) + Cm <i W + Cm q 2V + Cm 6 e ( 6e6e t )] 


q = p 


y 2 Sc l 

2Iy 


( 1 ) 

( 2 ) 

(3) 


for the lateral equations: 

2 / 

V = -ru + pw + g COS e sin 0 + p^.\C Yf0 + Cy^/3 + Cy r ||r + C Yp || + Cy^Sj 

P - * V + C 'r i + % & • * <V a + °* 6l >) 

i- = Pli^(°n,o + c ty/ + c ip SH"* c % fir* c ns a c a + c n 5 r «r) 


(4) 

(5) 

( 6 ) 


Because of the nature of the control inputs, it was possible to separate the longitudinal 
and lateral modes of motion; therefore, the parameter -extraction program solved for the 
longitudinal and lateral aerodynamic parameters separately. 


APPLICATION OF EXTRACTION PROGRAM 


Longitudinal Aerodynamic Parameters 

High-speed configuration .- Initial values of the various aerodynamic parameters 
which were required for the equations of motion were estimated by use of reference 6 
and are listed in table HE. The likelihood function initially included the states u, w, 
and q. An attempt was made to extract all the parameters listed in the longitudinal 
equations of motion and also the initial conditions for a best fit. (See run 1, table IV.) 
For this case, the mean-squared error decreased for several iterations and then began 
to increase. After about 15 iterations, the extracted values and the covariance matrix 
for the unknown parameters were read out. 
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Examination of the printout showed that at termination, changes in some of the 
unknown parameters were still large. (See run 1, table V.) The parameters with large 
variations from iteration to iteration changed in pairs or groups of three. Also, most of 
the parameters which were still varying had values significantly different from those pre- 
dicted by aerodynamic theory. Some groups of parameters were varying together and 
had high estimated correlation coefficients in the covariance matrix. (See table VI.) 
However, since the problem had not converged, the numbers in the covariance matrix 
could only be considered estimates of the actual covariances. In this problem, it was 
found that when the estimated correlation coefficients were near 1, the coefficients 
involved tended to deviate toward questionable values. The parameters which appeared 
to be related were C m q and C m( j; C mq and Cmg e J ^m a and C m ^. Examination 

of the flight data (fig. 5) showed that q, w, a z , and 5 e were all damped sine waves 
and approximately in phase so that a relation between these parameters might be 
expected. 

Since the pairs of coefficients which appeared to be related were coefficients which 
varied together from iteration to iteration, fixing one of each pair at a specified value 
could possibly force the other parameter in the pair to a specific value. Therefore, 

C m ^ was set equal to zero so that C m q assumed the whole effect of the C m q and 
C m< j combination and C m g was set equal to the value used in reference 1 which was 

taken from reference 9. Fixing the two coefficients reduced the number of unknowns 
being extracted. It was felt that if the number of unknowns were reduced further, there 
would be a better chance for convergence, and then more unknowns could be added. There 
fore, the forward velocity equation was deactivated, and the initial conditions of w and 
q were considered to be those of the flight data before any controls were activated, and 
Cz, 0 and C m>0 were calculated to balance the equations initially. Therefore, u Q , 

Cx,o Cx a > w o> C z>0 , q 0 , and C m>0 also were eliminated as unknowns in the mathe- 
matical model. (See run 2, table IV.) The program was rerun. The mean-squared error 
stabilized and became constant after about 10 iterations. 


By using these newly selected parameters as starting values, the initial conditions 
on u, w, and q, and the parameters Cx >0 > C Zo , C m>0 , and Cx^ were added as 
additional unknowns. (See run 3, table IV.) Again, this program stabilized to a lower 
performance index and a set of parameters was obtained. At this point Cx f o> ^z o> 
and C m>0 seemed to be well determined. A set of parameters had been obtained but 
c mg e and C m ^ had been set to specific numbers. Also, the value obtained for C Zg 

had a standard deviation of about 100 percent; the value obtained for was about 

85 percent higher than expected (table IH); and C mQ , was about one -third the value 
expected. 
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These results were not felt to be satisfactory. The large standard deviation in 

C z - indicated that this parameter could not be determined accurately from the avail - 
°e 

able data. It was, therefore, decided to use the geometric relationship 



( 7 ) 


to estimate 


and to keep that value constant during further parameter extractions. 


In addition, previous unpublished results of a derivative -extraction study had indicated 
that inclusion of a parameter Cz q in the mathematical model might have an appreciable 
effect on C Zq; . The parameter C z ^ therefore was added to the model, and equation (2) 

became (see also run 4, table IV) 


w = g cos 9 + qu + j P 


c Z,o + C Z, 


, 6 e “ 6 e t 


Adding Cz q to the model changed Cz a from -6.89 to -5.5, which appeared to be more 
reasonable, but probably is still too negative. The parameter Cz q was found to be 
-26.9; however, it had a large estimated standard deviation. (See run 2, table V.) In 
addition, the estimated standard deviation of w was high, 1.25 m/sec (4.1 ft/sec) 

(run 2, table V) compared with the flight range of from about 3.048 to -6.096 m/sec 
(10 to -20 ft/sec) (fig. 5). Most of this difference is attributable to the uncertainty in 
measured angle of attack. The uncertainty was ±0.8°, which resulted in an uncertainty 
of 0.975 m/sec (±3.2 ft/sec) in w. It appeared advisable, therefore, to expand the like- 
lihood function to include additional flight data. In this case, the measured normal 
accelerations appeared to be pertinent and reasonably accurate and were therefore added 
to the likelihood function. 

After az was put into the likelihood function (in place of w) along with u and 
q, the value of C Zo , changed from -5.5 to -4.33, the values of the other unknown param- 
eters c hanging less than 5 percent. At this point u, w, q, and a z were all used in 
the likelihood function, and the program was run again. (See run 3, table V.) The values 
previously obtained changed by less than 2 percent. The program weighted the a z data 
so that the u, q, a z likelihood function gave about the same results as the u, w, q, 
a z likelihood function. At this point the estimated derivatives appeared to be reason- 
able, except possibly for C mQ , which was somewhat less negative than had been expected. 
(Compare table HI with run 3 of table V.) 

It should be noted that for runs 2 and 3 of table V, the parameter C m ^ had been 
set equal to zero, with the expectation that the extracted C mq would reflect any C m ^ 
contribution to the aircraft motion. As a matter of interest, it was decided to now insert 
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and hold constant a value of C m ^ equal to the value given in table HI. The results from 
the extraction program are given in table VII. A comparison of run 3 of table V and the 
values of table VII show that 


(1) The value of C mQ! changed from -0.384 to -0.63, which is close to the value 
extracted in the study of reference 1. 

(2) The value of C m changed from -24.7 to -18.1. However, note that the sum 

H 

Cmq + C m ^ is about the same for runs 2 and 3 of table V and the values given in 
table VH. The other parameters were not affected. The standard deviations of the 
parameters determined, except ^x a ’ were l ess than 15 percent of the value of the 
parameter and indicated that they were well determined for the model being used. The 
fit obtained to the flight data is shown in figure 5 and the derivatives are given in table VH. 


As a matter of interest, it was decided to make additional runs to examine the sig- 
nificance of the high estimated correlations indicated in table VI. These additional runs 
held fixed values of one of the parameters involved and allowed the program to extract the 
remaining parameters. The results were then examined to note how the remaining 
parameters changed as a function of the fixed parameter. In addition, the changes in the 
determinant of R were also observed. Some of the results are shown in figures 6 and 7. 
The parameters shown in the figures are those that varied by more than about 5 percent 
as the "fixed" parameter was changed about ±100 percent. Figure 6 shows that since the 
determinant of R remained very nearly constant over the range of C m ^ values, and 


since the extracted C 


m 


a 


and 


approximately true linear correlation between C m and 

4 


Cmq values varied linearly with C m ^, there exists an 


C m ^ and between C mQ , and 


C m ^. However, as shown in figure 7, the determinant of R varied somewhat as C m6 
was changed over the range shown. Therefore, a linear correlation between C mq and 
Cmg e or between Cgq and C mge cannot be definitely established. 

The preceding discussion illustrates the considerations necessary to obtain coeffi- 
cients for the set of data used herein. These considerations are 


(1) Obtain as good a set of starting values as feasible 

(2) Even if the program has not truly converged, the estimated covariance matrix 
for the unknown parameters will show potentially related parameters. 

(3) If the program is slowly diverging because of pairs of parameters varying, 
fixing one of the parameters based on other factors such as wind-tunnel data or refer- 
ence 6 will generally cause the program to converge to a set of derivatives. 
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(4) Once the performance index has stabilized for a particular set of parameters, 
the states in the performance index and the parameters in the mathematical model can be 
changed to determine whether any improvements in convergence occurs. 

(5) After investigating the importance of coefficients which could affect the calcu- 
lated response, and determining a minimum performance index, curves to show possible 
relations between coefficients can be established. 

(6) Based on all this information, a set of parameters for the airframe investigated 
can be determined. 

Low-speed configuration.- The general procedures given for the high-speed config- 
uration were followed for the low-speed configurations. The final values of the extracted 
parameters are shown in table VII. Time histories of the measured and calculated 
motions are shown in figure 8. 

Lateral Aerodynamic Parameters 

High-speed configuration.- There were two sets of flight records from which the 
lateral aerodynamic parameters could be obtained: the response to an aileron doublet, 
and the response to a rudder pulse. As will be indicated, both sets of flight records were 
used in determining the parameters. 

The initially estimated lateral parameters required to start the computing and 
extraction process were computed by using reference 6 and are shown in table HI. The 
program gave a fair match to the aircraft motions in about six iterations; however, the 
mean- squared error, after decreasing for three iterations, began to slowly increase. 

The derivatives extracted for this case (high-speed, aileron doublet input) are indicated 
in table VHI (run 1). Since the standard deviations were small for all parameters, with 
the exception of Cyp and Cy r , the parameters appeared to be well defined. (See 
table EX.) Varying Cyp and Cy r by ±100 percent did not affect the values of the other 
parameters or the fit, and since they were not important, they were subsequently dropped 
from the mathematical model. 

The off-diagonal terms of the estimated covariance matrix are shown as table X. 

A possible relation appeared to exist between C? and and between C n and 

Cng • In a similar situation for the longitudinal case, one of the apparently related var- 
iables assumed to be known, resulted in a nonincreasing mean-squared error after sev- 
eral iterations. For the lateral runs, it was thought that the rudder pulse input data 
could be used to determine C n p and Cjp. However, when rudder data were run (run 2, 
table VTEI), the mean-squared error began to increase rapidly after one or two iterations. 
The printout showed that the roll and sideslip parameters had large variations and the 
estimated covariance matrix implied a possible relation between C n p and Cn^, and 
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between and C^. Because these possible relations prevented a nonincreasing 

m£an-squared error, a procedure for using both sets of lateral data was developed. The 
procedure is as follows: 


(1) Obtain the control parameters and C n g from wind-tunnel tests (ref. 9), 

and hold these values constant during the following parameter extraction process. Fix 
the initial conditions at the flight values and calculate the trim coefficients to initially 
zero the equations, and fix the trim coefficients at these values. (See run 3, table VTLL) 
Extract the remaining parameters. By using the initial conditions and these parameters 
as initial estimates, rerun the problem with the initial conditions and trim coefficients as 
additional unknowns. (See run 4, table VTEI.) 


(2) A similar procedure is used with the rudder input flight data. The initial condi- 
tions and trim were first assumed to be known, and C/„ and C n „ were fixed at the 
values obtained from run 4. With these conditions, the problem was run until a non- 
increasing mean-squared error was obtained. (See run 5, table VIII.) The extracted 
parameters were then used as initial estimates, and the initial conditions and the trim 
coefficients were assumed to be unknown and the run repeated. (See run 6, table VIII.) 

(3) The results of runs 4 and 6 were compared and it was found that the differences 
in the parameters obtained varied less than 10 percent with the exception of Cn r which 
had a variation slightly greater than 10 percent. 

(4) The aileron input case was run with the yawing and sideslip derivatives obtained 
from the rudder input runs. (See run 7, table VIE.) The rudder input case was then run 
with the roll, sideslip, and yaw derivatives obtained from the aileron input runs. (See 
run 8, table VIII.) Runs 7 and 8 were repeated several times and a set of derivatives 
which best fit both sets of data was determined. These values are given in table XI and 
the resulting motion time histories are given in figures 9 to 12. 


As was done in the longitudinal runs, curves showing any significant percent changes 
in the parameters or in the determinant of R due to changes in the parameter assumed 
to be known are shown as figures 13 to 15. As before, in addition to the parameters pri- 
marily affected, there was some effect on other parameters. Therefore, the determi- 
nant of R did not remain constant and correlations could not be determined. 


DISCUSSION OF RESULTS 


The results shown in figures 5, 8, 9, 10, 11, and 12 show that the parameter- 
extraction program was able to obtain a set of aerodynamic parameters which permitted 
estimates of the aircraft motion compatible with the accuracies of the measurements of 
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the flight data. The results of this study are compared with results of references 1, 8, 
and 9 in tables XII and XHI for the longitudinal and lateral parameters, respectively. 
Major differences in values of the parameters are discussed. In this discussion the Tg 
for the wind tunnel was approximately the same as the Tg for flight. 

Longitudinal Parameters 

The longitudinal parameters obtained from the present study and those of refer- 
ences 1, 8, and 9 are in reasonable agreement in most instances. However, there are 
several important differences which require discussion. 

Cx a -~ There was not enough variation in forward velocity to permit accurate deter- 
mination of C Xq! from the flight data. 

All references had about the same value of exce Pt that for reference 1 

which is somewhat higher than appears to be reasonable. In searching for reasons for 
this difference, it was noted that the mathematical model used in reference 1 did not 
include the parameter C z , whereas it was included in the present study. The param- 

VI 

eter C z , therefore, was set equal to zero in the present study, and the program was 

4 

used to extract the remaining parameters. It was found that most of the parameters 
remained at the values of table XII; however, C Za increased (negatively) from -4.33 
to about -6.0. It appears that, at least for the Navion, it was important to include C z 

VL 

among the parameters in the mathematical model. 

Cm tf - The present study yielded a value of C mQ , which was less than one -half the 

wind-tunnel value (ref. 9) if C m ^ were set to zero in the mathematical model. How- 
ever, if C m ^ were set equal to -6.5 for the high-speed case and to -6.0 for the low- 
speed case, C mQ! was about 25 percent less than the value obtained from wind-tunnel 
tests (ref. 9). The difference between the wind-tunnel value of C mQ! and the value 
extracted in the present study raised a question of uniqueness of derivatives. It was, 
therefore, decided to insert the wind-tunnel values of the parameters as constants in the 
parameter -extraction program, and to let the program extract the remaining parameters. 
The results of the time-history calculation are shown in figures 16 and 17, and the fixed 
wind-tunnel parameters and the remaining (extracted) parameters are listed in table XIV 
for the high-speed and low-speed configurations, respectively. Although the computed 
time histories are in fair agreement with the flight data, the agreement is not as good as 
that shown in figures 5 and 8. Also note that the period of the calculated time histories 
is somewhat shorter than that of the flight data (figs. 16 and 17), which is a consequence 
of increasing the value of C mQ , to the wind-tunnel value. 
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Lateral Parameters 





For the high-speed lateral runs, the agreement between the results of reference 1 
and the present study was generally good. (See table XIII.) The major exception was a 
difference in C n ^ of about 15 percent. Also C^ r showed a 40-percent difference and 
C ng a 20-percent difference. However, the value obtained for C^ r in reference 1 was 

based on flight data different from those used in this study. More serious differences 
arose between the present study and reference 9 where Cy^,, C^, and C n ^ were all 
different by about 40 percent. Significant differences exist between the values of Cy^, 

C and C ng ^ obtained in reference 9 and those of the present study. 

By examining the low -speed lateral runs, differences between the results of refer- 
ence 1 and the present study were found in C ir (100 percent), C nr (35 percent), and 
Cng (30 percent). Differences between reference 9 and the present study occurred in 
C ip (90 percent), C n ^g (50 percent), and C^. 

CONCLUDING REMARKS 

An iterative procedure, which is characterized as a maximum-likelihood minimum- 
variance technique, was used to extract the aerodynamic parameters of a Navion airplane 
from flight data. Results from the study showed that the parameter -extraction program 
can produce aerodynamic parameters which will permit close estimation of the aircraft 
time histories used in the extraction process. The study also showed that good judgment 
and ingenuity is required to circumvent problems occurring because of apparent relations 
between some of the parameters. It also became evident that the mathematical model 
selected to represent the aircraft is very important. 

The extracted parameters agreed reasonably well with the values obtained by analog 
matching the same data, with the exception of the change in normal -force coefficient with 
angle of attack The agreement with wind-tunnel parameters was not as good for 

the variations of pitching-moment coefficient with angle of attack side-force coef- 

ficient with sideslip angle (Cy^, rolling-moment coefficient with sideslip angle (C^, 
and yawing-moment coefficient with sideslip angle However, of the parameters 

determined by the program, only one had a standard deviation greater than 15 percent of 
the value of the parameter and the parameters gave a reasonable fit to the flight data. 

The parameter estimation process was greatly aided by the capability of the opera- 
tor to interact with the program. The hands-on feature was especially useful during the 
initial stages of the estimation procedure when a number of different parameters were 
being considered as known or unknown in the mathematical model. Some of the other 
advantages of this parameter estimation program were found to be: 
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(1) The mathematical logic for updating unknown parameters is mechanized. 

(2) A "best fit," within the limits of the mathematical model and the data used, is 
obtained. 

(3) The estimated standard deviation of each unknown parameter is obtained. 

The control inputs used to obtain the aircraft responses used herein caused these 
responses to be approximately sine waves which were close to being in phase. Uniquely 
determining the parameters using this flight data was difficult and the values of several 
parameters had to be assumed to be known from theoretical calculations or wind-tunnel 
data. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., January 20, 1972. 
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TABLE I.- NAVION AIRPLANE DIMENSIONS 


Wing; 

Area, S, m2 (ft2) 17.112 (184) 

Sweep, leading edge, deg 2.996 

Aspect ratio, A 6.04 

Taper ratio, X 0.54 

Mean aerodynamic chord, m (ft) 1.74 (5.7) 

Dihedral, deg 7.5 

Incidence at root, deg +2 

Incidence at tip, deg -1 

Airfoil: 

Tip NACA 6410 R 

Root NACA 4415 R 

Horizontal tail: 

Area, m2 (ft 2 ) 4.0 (43) 

Sweep, leading edge, deg 6 

Aspect ratio 4.0 

Taper ratio 0.67 

Airfoil NACA 0012 

Incidence, deg -3 

Vertical tail: 

Area (above horizontal stabilizer), m2 (ft 2 ) 1.163 (12.5) 

Airfoil: 

Root Modified NACA 0013.2 

Tip Modified NACA 0012.04 

Fin offset, deg 2 

Propeller characteristics: 

Diameter, m (in.) 2.14 (84) 

Number of blades 2 

Side-force factor 100 

Power plant: 

Continental engine Model no. 10520B 

Horsepower rating at take-off at 2700 rpm 285 

Mass and inertia characteristics for data of this report: 

Gross mass, kg (lb) 1335.76 (2948) 

Center of gravity, percent c 25 

I X , kg-m 2 (Slug-ft 2 ) 1742.33 (1284.08) 

I Y , kg-m 2 (slug-ft 2 ) 3762.4 (2772.86) 

I z ,kg-m 2 (slug-ft 2 ) 4389.10 (3234.72) 

Control surfaces: 


Control surface 

Area 

Deflection, deg 

m 2 

ft 2 

Flaps (plain) 

7.775 

83.6 

40 

Stabilizer 

2.79 

30.0 


Elevator 

1.31 

14.1 

Up 30, down 20 

Aileron 

.502 

5.4 

20 

Rudder 

.558 

6.0 
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TABLE H.- FLIGHT-TEST CONDITIONS FOR THE VARIABLE-STABILITY NAVION 


Condition I 

Altitude, km (ft) 1.525 (5000) 

Velocity (true airspeed), m/sec (ft/sec) 73.2 (240) 

Flaps, deg 0 

Lift coefficient 0.271 

Thrust coefficient (estimated), Tj, 0.02255 


Condition H 

1.525 (5000) 
43.9 (144.1) 
20 
0.753 
0.05716 
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TABLE ID.- INITIAL ESTIMATES OF THE 


AERODYNAMIC COEFFICIENTS* 



Value 

0.0015 

0 

0.02 

- 5.23 

- 0.40 

0 

- 1.26 

- 6.50 

- 13.00 

- 1.4 

- 0.35 

0.50 

- 0.30 

0.08 

- 0.06 

0.04 

- 0.45 

0.15 

0.30 

0.049 

- 0.024 

- 0.082 

- 0.005 

- 0.08 


^Estimated by using methods presented 
in reference 6. 





TABLE V.- PARAMETERS DETERMINED AT VARIOUS STAGES OF FITTING 
THE LONGITUDINAL FLIGHT DATA 









Run after adding a^ 

to the 





Hun after establishing better 

oerformance index to the 



Hun 1 

starting values and adding 
Cry to the model (Run 2) 

model and fixing C 

m «e 

Parameter 




q 



and 

C z (Run 

u e 

3) 


Value 

AP 

Standard 

deviation 

Value 

AP 

Standard 

deviation 

Value 

AP 

Standard 

deviation 

Cv 

A a 

-0.69 

-0.15 

0.04 

-0.09 

6.6 X 10" 6 

0.04 

-0.265 

5.3 X 10-6 

0.05 

C7 

^ a 

-6.87 

0.16 

0.12 

-5.47 

6.0 X 10- 4 

0.15 

-4.33 

7 x 10-5 

0.007 

° Z q 

0 

0 

— 

-26.9 

. . 

0.1 

19.6 

-15.9 

-4.23 x IO- 4 

2.28 

Cz 5e * 

-4.11 

0.33 

0.15 

-0.511 

0 

— 

-.0.511 



Cm a 

3.00 

-0.05 

0.43 

-0.584 

2.3 x 10- 7 

0.02 

-0.39 

2 x 10- 7 

0.011 


124.9 

-1.2 

15.5 

0 

0 


0 



Cm a 





Cm q 

-144.2 

1.18 

16.3 

-24.7 

9 X 10" 5 

0.43 

-24.5 

8 x 10-5 

... 

0.41 


-0.91 

-0.23 

0.14 

1.42 

0 


-1.42 



C m 6 

°e 




Standard deviation 










of u 

28 




4.4 


1.7_ 



Standard deviation 






1 




of w 

8 




4.1 


2.5 



Standard deviation 










of q 

0.003 



1 

0.0004 


L3,7_x_ig- 5 



Standard deviation 










of a z 







1.6 X 10-3 





24 



TABLE VI.- ESTIMATED CORRELATION MATRIX FOR ELEVATOR DOUBLET 

INPUT FOR RUN 1 OF TABLE V 

[V = 73.2 m/sec (240 ft/sec)J 



Cx 

A a 

C Z 

Zq, 

Cz 6 e 

c m 

a 

c m • 

a 

o 

3 

►a 

Cm 6 e 

CXa 

1 

0.08 

0.15 

-0.04 

0.03 

0.015 

-0.13 

c z a 

0.08 

1 

0.69 

-0.11 

0.09 

-0.06 

-0.10 

C Z6 e 

0.15 

0.69 

1 

0.07 

0.10 

-0.11 

0.034 

c m a 

0.04 

-0.11 

0.07 

1 

0.978 

0.762 

0.15 

° m a 

0.03 

0.09 

0.10 

0.978 

1 

0.997 

-0.014 

c m q 

0.015 

-0.06 

-0.11 

0.762 

0.997 

1 

0.968 

Cm 6 e 

-0.13 

-0.10 

0.034 

0.15 

-0.014 

0.968 

1 
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TABLE VH.- LONGITUDINAL AERODYNAMIC PARAMETERS OF THE 
NAVION AIRPLANE EXTRACTED FROM FLIGHT DATA 


(a) 

V Q = 73.2 m/sec (240 ft/sec); 5 f = 0° 


Parameter 

Value 

Standard deviation 



0.262 

0.046 

C Z a 

-4.33 

0.007 

C Z q 

-15.9 

2.28 

c z 6e 

-0.511 


c m a 

-0.63 (-0.77 corrected to 25% c.g.) 

0.006 

Cm ci 

-6.5 

0.52 

Cm q 

-18.1 


Cm 6 e 

-1.42 


(b) 

V D = 43.9 m/sec (144.1 ft/sec); 6f = 20 

o 

Parameter 

Value 

Standard deviation 

c x« 

1.37 

0.09 

O 

N 

Q 

-4.86 

0.09 

C Z q 

-27.13 

1.75 

C Z 6e 

-0.52 


Cma 

-0.70 (-0.84 corrected to 25% c.g.) 

0.011 

Cm a- 

-6.0 

0.25 

O 

3 

-16.4 

0.035 

Cm 6 e 

-1.55 
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TABLE Vffl.- PARAMETERS USED IN MATHEMATICAL MODEL FOR 

LATERAL RUNS 

[Key: - denotes parameter assumed known; y parameter assumed 
unknown; * not in mathematical model. Subscript 1 indicates 
aileron-input data; subscript 2 indicates rudder -input data] 

Parameter Run lj I Run 22 Run Run 4j_ Run 52 Run 62 Run 7^ I Run 82 


c Y,o 

% 

C Yp 

C Y r 

C Y 6r 












TABLE IX.- AERODYNAMIC PARAMETERS FOR NAVION AIRPLANE 
EXTRACTED FROM AILERON DOUBLET RUN 

[v = 73.2 m/sec (240 ft/sec); 5 f = 0°] 


Parameter 


Value 


Standard deviation 



TABLE X.- ESTIMATED CORRELATION MATRIX FOR AILERON DOUBLET FOR 
V = 73.2 m/sec (240 ft/sec) AND 6f = 0 ° 

^M^R^m) with the diagonal terms set to lj 



Cy 

Y /3 

Cy 

*P 

c v 

x r 

% 

Ci 

l P 

C, 

l r 


Cn /3 

c n p 

c n r 

Cn 6 

°a 

C Y/3 

1 

- 0.17 

- 0.14 

0.02 

- 0.003 

0.31 

- 0.03 

0.03 

- 0.04 

- 0.84 

- 0.10 

°Yp 

- 0*17 

1 

0.23 

0.10 

0.33 

0.004 

- 0.30 

- 0.28 

- 0.61 

0.13 

- 0.65 

C Y r 

- 0.14 

0.23 

1 

- 0.34 

0.13 

0.001 

- 0.11 

0.76 

- 0.22 

0.11 

- 0.22 

% 

0.02 

0.10 

- 0.34 

1 

0.18 

0.04 

- 0.18 

- 0.07 

- 0.02 

0.04 

0.04 

% 

- 0.003 

0.33 

0.13 

0.18 

1 

0.47 

- 0.94 

- 0.02 

0.12 

0.12 

0.08 

C lr 

0.31 

0.004 

0.001 

0.04 

0.47 

1 

- 0.35 

- 0.07 

0.19 

0.11 

0.10 


- 0.03 

- 0.30 

- 0.11 

- 0.18 

- 0.94 

- 0.35 

1 

0.02 

- 0.02 

- 0.03 

0.03 

c np 

0.03 

- 0.28 

0.76 

- 0.07 

- 0.02 

- 0.07 

0.02 

1 

0.14 

- 0.01 

0.15 

c n p 

- 0.04 

- 0.61 

- 0.22 

- 0.02 

0.12 

0.19 

- 0.02 

0.14 

1 

0.26 

0.97 

Cnr 

- 0.84 

0.13 

0.11 

- 0.04 

0.12 

0.11 

- 0.03 

- 0.01 

0.26 

1 

0.25 

Cn 6a 

- 0.10 

- 0.65 

- 0.22 

0.04 

0.08 

0.10 

- 0.03 

0.15 

0.97 

0.25 

1 
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TABLE XI.- LATERAL PARAMETERS 


Parameter 

V = 73.2 

m/sec (240 ft/sec); 

6f = 0° 

V= 43.9 

m/sec (144.1 ft/sec); 
6 f = 20° 


Value 

Standard deviation 

Value 

Standard deviation 

Cy 

Y /3 

-0.6 

0.078 

-0.74 

0.059 

°Ya r 

0.33 

0.05 

0.68 

0.03 


-0.07 

0.0016 

-0.053 

0.0012 

C h 

-0.49 

0.003 

-0.53 

0.0037 

C h 

0.11 

0.008 

0.114 

0.0068 

Cl 

5 a 

0.154 

0.0004 

0.16 

0.0004 

C7 

0.026 

0.0018 

0.0007 

0.0013 

Cn /3 

0.073 

0.0004 

0.083 

0.0005 

c n p 

-0.04 

0.0005 

-0.147 

0.0016 

Cn r 

-0.09 

0.005 

-0.108 

0.0044 

c % 

-0.004 

8 x 10-6 

-0.0015 

0.0001 

Cn 6 r 

-0.063 

0.0007 

-0.067 

0.0005 
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TABLE XH.- COMPARISON OF LONGITUDINAL PARAMETERS 
WITH THOSE OF REFERENCES 


Parameter 

Present study 

Reference 1 

Reference 9 

Reference 8 



6 f = 0° 



C Xa 

0.262 


0.12 

0.33 

c z a 

-4.33 

-6.04 

-4.52 

-4.44 

C Z q 

-15.9 




C ^a 

-0.511 


-0.43 

-0.355 

C m a 

-0.77 


-0.95 


Cm q + 

-24.6 

-18.3 


-14.32 

C “5 e 

-1.42 

-1.42 

-1.42 




6 f = 20° 



C x a 

1.37 


0.30 


o 

N 

ft 

-4.86 

-6.40 

-4.52 


° Z q 

-27.13 




Cz 6e 

-0.52 


-0.53 


C % 

-0.84 




C m q + c m,j 

-22.4 

-15.5 



Cm «e 

-1.55 

-1.55 

-1.55 

i 




i i ■!■■■■ us min ■■ 


TABLE m- COMPARISON OF LATERAL PARAMETERS WITH THOSE OF REFERENCES 
Present study 



- 0.6 

0.33 

-0.07 

-0.49 

0.11 

0.154 

0.026 

0.073 

-0.04 

-0.09 

-0.004 

-0.063 

-0.74 

0.68 

-0.053 

-0.53 

0.114 

0.16 

0.007 

0.080 

-0.147 

- 0.12 

-0.0015 

-0.075 


Reference 1 
6 f = 0° 
-0.61 


-0.067 

-0.46 

0.07 

0.152 

0.086 

-0.038 

-0.088 

-0.0047 

-0.075 
6 f = 20° 


-0.051 

-0.48 

0.27 

0.150 

0.084 

-0.141 

-0.163 

-0.0013 

-0.093 


Reference 9 

- 0.8 

0.143 

-0.109 


0.152 

0.029 

0.109 


-0.0047 

-0.077 


-0.77 

0.143 

-0.096 


0.09 

0.015 

0.126 


Reference 8 

-0.56 

0.157 

-0.074 

-0.41 

0.107 

0.134 

0.012 

0.070 

-0.058 

-0.125 

-0.0035 

-0.07 


0 

-0.08 


i 

-i 

£ 
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TABLE XIV.- VALUES OBTAINED BY FIXING PARAMETERS OBTAINED 


FROM REFERENCE 9 AND DETERMINING C mq + C mA 


(a) V Q = 73.2 m/sec (240 ft/sec); 6f = 0° 



Parameter 


Value 

Standard deviation 

C y . . . • • 


• •it* 

0.12 


Cz* 


• • • • • 

-4.52 


° Z q 


• • • • • 

-5.25 


Cz 6e 


• • • • « 

-0.403 


C % 


• • • • • 

-0.95 


Cmq + Cm^ • 


• • • • • 

-15.6 

1.63 

c m g 

°e 



-1.42 



(b) V G = 43.9 m/sec (144.1 ft/sec); = 20° 


r~ 



Parameter 

Value 

Standard deviation 

c x<* 

0.3 


c z a 

-4.52 


C Zq 

-6.06 


Cz 6 e 

-0.53 


C ma 

-0.95 


C m q + C m(i 

-13.12 

1.31 

Cm ®e 

-1.55 
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Figure 2.- Three -view drawing and principal dimensions. All dimensions are in feet 
(meters). Tail incidence angle ij- is normally set at 0°; it was set at -5° for these 
tests. 
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Figure 3.- Flow diagram for parameter -extraction program. 




Flight and computed data 





Time, sec 

Figure 5.- Comparison of flight data with time histories computed by using the aerody- 
namic parameters of table VII for elevator doublet input. V = 73.2 m/sec (240 ft/sec); 
6 f = 0°. Flight and computed data are the same for control inputs. 
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Figure 7.- Percent changes in the values of C m q, C mo! , and the determinant of R 

for changes in the value of Cm* . 
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Figure 8.- Comparison of flight data with time histories computed by 
using aerodynamic parameters of table VII for elevator doublet input. 
V = 43.9 m/sec (144.1 ft/sec); = 20°. Flight and computed data 
are the same for control inputs. 
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Figure 9.- Comparison of flight data with time histories computed by 
using aerodynamic parameters of table XI for aileron doublet input. 
V = 73.2 m/sec (240 ft/sec); 6 f = 0°. Flight and computed data 
are the same for control inputs. 
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Figure 10.- Comparison of flight data with time histories computed by 
using aerodynamic parameters of table XI for rudder pulse input. 

V » 73.2 m/sec (240 ft/sec); = 0°. Flight and computed data 
are the same for control inputs. 













radians 8 e , radians 8 n , radians 


Flight and. computed data 


Time, sec 

Figure 11.- Comparison of flight data with time histories computed by 
using aerodynamic parameters of table XI for an aileron doublet input 
V = 43.9 m/sec (144.1 ft/sec); 6 f = 20°. Flight and computed data 
are the same for control inputs. 
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Figure 11.- Concluded. 
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Figure 12.- Comparison of flight data with time histories computed by 
using aerodynamic parameters of table XI for a rudder pulse input. 
V = 43.9 m/sec (144.1 ft/sec); 6 f = 20°. Flight and computed 
data are the same for control inputs. 
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Figure 12.- Concluded. 




Figure 13.- Percent changes in the values of C^, C , and the determinant of R 

for changes in the value of C^ p . 








Figure 14.- Percent changes in the values of C np , C lr , C^, and the determinant of R 

for changes in the value of Cn 5 . 
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Figure 16.- Response of the mathematical model using derivatives from 

reference 9 and from the extraction program as compared with flight data 
for the 0° flap, V = 73.2 m/sec (240 ft/sec) case when the control input 
was an elevator doublet. 
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Figure 17.- Response of the mathematical model using derivatives from 
reference 9 and from the extraction program as compared with flight 
data for the 20° flap, V = 43.9 m/sec (144.1 ft/sec) case where the 
control input was an elevator doublet. 
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Figure 17.- Concluded. 
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